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We present a detection of the unnormalized skewness (T^(n)) induced by the thermal Sunyaev- 
Zel'dovich (tSZ) effect in filtered Atacama Cosmology Telescope (ACT) 148 GHz cosmic microwave 
background temperature maps. Contamination due to infrared and radio sources is minimized by 
template subtraction of resolved sources and by constructing a mask using outlying values in the 
218 GHz (tSZ-null) ACT maps. We measure (T^(n)) = -31 ±6 /iK^ (Gaussian statistics assumed) 
or ±14 /xK^ (including non-Gaussian corrections) in the filtered ACT data, a 5cr detection. We show 
that the skewness is a sensitive probe of as, and use analytic calculations and tSZ simulations to 
obtain cosmological constraints from this measurement. From this signal alone we infer a value of 
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0.79 
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(68% C.L. 
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-0.06 



(95% C.L.). Our results demonstrate that measurements of non- 



Gaussianity can be a useful method for characterizing the tSZ effect and extracting the underlying 
cosmological information. 
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I. INTRODUCTION 

Current observations of the cosmic microwave back- 
ground (CMB) anisotropics on arcminute scales using 
experiments such as the Atacama Cosmology Telescope 
[ACT; H-g and the South Pole Telescope [SPT; U, 15| 
probe not only the primordial microwave background 
fluctuations sourced 13.7 billion years ago, but also 
measure secondary anisotropics caused by more recent 
and less distant physical processes. Such secondary 
anisotropics are induced by infrared (IR) dusty galax- 
ies and radio sources, gravitational lensing, and the 
Sunyaev-Zel'dovich (SZ) effect. The SZ effect %^ arises 
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due to the inverse Compton scattering of CMB photons 
off high energy electrons located predominantly in hot 
gas in galaxy clusters (the intra-cluster medium, or ICM). 
This scattering modifies the spectrum of CMB photons in 
the direction of a cluster in a way that depends on both 
the thermal energy contained in the ICM (the thermal 
SZ effect) as well as the peculiar velocity of the cluster 
with respect to the CMB rest frame (the kinetic SZ ef- 
fect). The kinetic SZ effect simply increases or decreases 
the amplitude of the CMB spectrum in the direction of 
a cluster, but the thermal SZ (tSZ) effect modifies the 
CMB spectrum in a frequency-dependent manner. The 
tSZ effect is characterized by a decrease (increase) in the 
observed CMB temperature at frequencies below (above) 
218 GHz in the direction of a galaxy cluster due to in- 
verse Compton scattering. The thermal effect is generally 



at least an order of magnitude larger than the kinetic 
effect for a typical massive cluster at 148 GHz. Mea- 
surements of the tSZ signal, which is proportional to the 
integrated ICM pressure along the line of sight, can be 
used to observe the high redshift universe, constrain cos- 
mo logical parameters - in particular ag, the variance of 
matter fluctuations on scales of 8 Mpc/h - and probe 
baryonic physics in the ICM. 

The tSZ signal has so far primarily been studied either 
by directly resolvin g in dividual clusters in arcminute- 
scale CMB maps |8l-[l2| or by measuring it statistically 
through its presence in the small-scale CMB power spec- 
trum [13, 14]. However, it is by no means obvious that 
the power spectrum is the best way to characterize the 
statistical properties of the tSZ field. Indeed, measuring 
the tSZ signal in the power spectrum is challenging be- 
cause there are many other sources of CMB power on ar- 
cminute scales: primordial CMB fluctuations, CMB lens- 
ing, instrumental noise, dusty star-forming IR galaxies, 
and radio sources. In order to disentangle these contri- 
butions to the power spectrum and isolate the amplitude 
of the tSZ signal, a sophisticated multifrequency analysis 
is required, which involves modeling the power spectrum 
contribution of each of these components in at least two 
frequency bands. 

In this paper we instead measure the tSZ signal using 
the unnormalized skewness of the filtered temperature 
fluctuation (T^(n)). This quantity has the significant 
advantage that, unlike measurements of the tSZ effect 
through the power spectrum, its measurement does not 
require the subtraction of Gaussian contributions, be- 
cause it is only sensitive to non- Gaussian signals with 
non-zero skewness. The primordial CMB (which is as- 
sumed to be Gaussian on these scales) and instrumen- 
tal noise (which is Gaussian) hence do not contribute 
to it. In addition, CMB lensing and the kinetic SZ ef- 
fect do not induce skewness (as they are equally likely 
to produce positive and negative fluctuations), and so 
do not contribute either. The primary contributions 
to this quantity are thus only the tSZ effect and point 
sources. These signals have a different frequency depen- 
dence. Furthermore, the tSZ signal contributes nega- 
tive skewness, whereas radio and IR point sources con- 
tribute positive skewness. These characteristics allow the 
tSZ signal to be effectively isolated and studied, as first 
pointed out in [15|]. 

Measurements of the skewness also possess significant 
advantages from an astrophysical perspective. A con- 
sistent problem plaguing studies of the tSZ power spec- 
trum has been theoretical uncertainty in the ICM elec- 
tron pressure profile [l6l-[l8j, especially in the low-mass, 
high-redshift groups and clusters that contribute much 
of the signal. As discussed in the following section in de- 
tail, the tSZ skewness signal is dominated by character- 
istically higher-mass, lower-redshift clusters than those 
that source the power spectrum signal. The ICM astro- 
physics for these objects is better constrained by X-ray 
observations and they are less sensitive to energy input 



from non-gravitational sources [Isl, [IqI- Thus, the theo- 
retical systematic uncertainty in modeling the tSZ skew- 
ness is correspondingly lower as well. In addition, at 148 
GHz, dusty star- forming galaxies are less prevalent in 
massive, low-redshift clusters (which contribute more to 
the skewness) than in high-redshift groups and clusters 
(which contribute more to the tSZ power spectrum) [2Q| . 
Thus, we expect the correlation between tSZ signal and 
dusty galaxy emission, which can complicate analyses of 
the tSZ effect, to be smaller for a measurement of the 
skewness. 

Moreover, the tSZ skewness scales with a higher power 
of (78 than the tSZ power spectrum amplitude. This re- 
sult is precisely what one would expect if the signal were 
dominated by higher-mass, rarer objects, as the high- 
mass tail of the mass function is particularly sensitive 
to a change in ag. This provides the prospect of tight 
constraints on cosmological parameters from the skew- 
ness that are competitive with constraints from the power 
spectrum. 

In this paper, we first explain the usefulness of the 
skewness as a cosmological probe by theoretically deriv- 
ing its scaling with as as well as the characteristic masses 
of the objects sourcing the signal. Subsequent sections of 
the paper describe how we measured this skewness in the 
ACT data. We describe how the ACT temperature maps 
are processed in order to make a reliable measurement of 
the unnormalized skewness due to the tSZ effect, and dis- 
cuss how contamination from IR dusty galaxies and radio 
point sources is minimized. We report the measurement 
results and discuss how the errors are calculated. Finally, 
we discuss the cosmological constraints and associated 
uncertainties derived from this measurement. 

We assume a flat ACDM cosmology throughout, with 
parameters set to their WMAP5 values ^] unless other- 
wise specified. All masses are quoted in units of M^/h^ 
where h = i:fo/(100kms-^ Mpc"^) and Hq is the Hubble 
parameter today. 



II. SKEWNESS OF THE TSZ EFFECT 

In this section, we investigate the N^^ moments of the 
pixel probability density function, (T^) = (T(n)^), fo- 
cusing on the specific case of the unnormalized skewness 
(r^y We show that the unnormalized skewness (T^) has 
a steeper scaling with as than the power spectrum ampli- 
tude and is dominated by characteristically higher-mass, 
lower-redshift clusters, for which the ICM astrophysics 
is better constrained and modeled. As explained earlier, 
these characteristics make tSZ skewness measurements a 
useful cosmological probe. 

In order to calculate the N^^ moment of the tSZ field, 
we assume the distribution of clusters on the sky can be 
adequately described by a Poisson distribution (and that 
contributions due to clustering and overlapping sources 



are negligible \22\). The N^^ moment is then given by 

(1) 

where T{0; M, z) is the tSZ temperature decrement at 
position on the sky with respect to the center of a 
cluster of mass M at redshift z: 



T{e-M,z) = g(v)TcMB 



(7t 



X j P,Ljp + d\{z)\e\^;M,z\dl, (2) 

where g{v) is the spectral function of the tSZ effect, 
dA{z) is the angular diameter distance to redshift z, and 
the integral is taken over the electron pressure profile 
Pe{r] M, z) along the line of sight. 

For a given cosmology, Eqs. ([T]) and (|2]) show that there 
are two ingredients needed to calculate the N^^ tSZ mo- 
ment (in addition to the comoving volume per steradian 
dV/dz^ which can be calculated easily): (1) the halo mass 
function dn{M, z)/dM and (2) the electron pressure pro- 
file Peiv] M, z) for halos of mass M at redshift z. We use 
the halo mass function of [23] with the redshift-dependent 
parameters given in their Eqs. (5)-(8). While uncertain- 
ties in tSZ calculations due to the mass function are often 
neglected, they may be more important for the skew- 
ness than the power spectrum, as the skewness is more 
sensitive to the high-mass exponential tail of the mass 
function. We estimate the uncertainty arising from the 
mass function by performing alternate calculations with 
the mass function of [24], which predicts more massive 
clusters at low redshift than [23] for the same cosmol- 
ogy. As an example, the predicted skewness calculated 
using the pressure profile of [16| with the mass function 
of [23] is ^ 35% lower than the equivalent result using 
the mass function of [24]. However, the derived scalings 
of the variance and skewness with ag computed using [24] 
are identical to those found below using [23|. Thus, the 
scalings calculated below are robust to uncertainties in 
the mass function, and we use them later to interpret 
our skewness measurement. However, we rely on cosmo- 
logical simulations to obtain predicted values of the tSZ 
skewness. We do not consider alternate mass functions 
any further in our analytic calculations. 

We consider three different pressure profiles from [16|, 
llTI . l25| in order to evaluate the theoretical uncertainty in 
the scaling of the tSZ skewness with ag. These profiles 
differ in how they are derived and in the ICM physics 
they assume. They thus provide a measure of the scatter 
in the scalings of the variance and skewness with ag due 
to uncertainties in the gas physics. 

Finally, in order to make a faithful comparison be- 
tween the theory and data, we convolve Eq. (J2]) with 
the Fourier-space filter described in the subsequent data 
analysis sections of this paper. In addition, we account 
for the 12(7 pixel fluctuation cutoff used in the data anal- 
ysis (see below) by placing each "cluster" of mass M 



and redshift z in the integrals of Eq. ([T]) in an ideal- 
ized ACT pixel and computing the observed temperature 
decrement, accounting carefully for geometric effects that 
can arise depending on the alignment of the cluster and 
pixel centers. If the calculated temperature decrement 
exceeds the 12<j cutoff, then we do not include this clus- 
ter in the integrals. These steps cannot be neglected, as 
the filter and cutoff reduce the predicted tSZ skewness 
amplitude by up to 95% compared to the pure theoreti- 
cal value [26]. Most of this reduction is due to the filter, 
which modestly suppresses the temperature decrement 
profile of a typical cluster; this suppression strongly af- 
fects the skewness because it is a cubic statistic. 

The analytic theory described above determines the 
scaling of the N^^ tSZ moment with ag. In particular, 
we compute Eq. ^ with A^ = 2, A^ = 3, and TV = 6 for 
each of the chosen pressure profiles while varying cr^. The 
scalings of the variance (A" = 2), the skewness (A" = 3), 
and the sixth moment (A" = 6, which we require for error 
calculation) with as are well-described by power laws for 
each of these profiles: (T^'^.e^ ^ erg'''. For the profile 
of jl6|, we find 0^2 = 7.8, aa = 11.1, and a^ = 16.7; 
for the profile of [13, we find 0^2 = 8-0, 0^3 = 11.2, and 
aQ = 15.9; and for the profile of [25|, we find 0^2 = 7.6, 
as = 10.7, and as = 18.0. Note that the scaling of the 
variance matches the scaling of the tSZ power spectrum 
amplitude that has been found by a number of other 
studies, as expected (e.g., [25|, IS^). The scaling of the 
unnormalized skewness is similar to that found by [28], 
who obtained 0^3 = 10.25. Also, note that the skewness 
scaling is modified slightly from its pure theoretical value 
[26] due to the Fourier-space filter and pixel fluctuation 
cutoff mentioned above. The overall conclusion is that 
the skewness scales with a higher power of as than the 
variance (or power spectrum). We use this scaling to 
derive a constraint on as from our measurement of the 
skewness below. 

In addition, we compare the characteristic mass scale 
responsible for the tSZ skewness and tSZ variance (or 
power spectrum) signals. Analytic calculations show 
that the tSZ power spectrum amplitude typically receives 
^ 50% of its value from halos with M < 2-3x10^"^ M0//1, 
while the tSZ skewness receives only ^ 20% of its am- 
plitude from these less massive objects. This indicates 
that the clusters responsible for the tSZ skewness signal 
are better theoretically modeled than those responsible 
for much of the tSZ power spectrum, both because mas- 
sive clusters have been observed more thoroughly, and 
because more massive clusters are dominated by gravita- 
tional heating and are less sensitive to non-linear energy 
input from active galactic nuclei, turbulence, and other 
mechanisms [18, 19]. We verify this claim when inter- 
preting the skewness measurement below, finding that 
the systematic theoretical uncertainty (as derived from 
simulations) is slightly smaller than the statistical error 
from the measurement, though still non- negligible. 



III. MAP PROCESSING 

A. Filtering the Maps 

The Atacama Cosmology Telescope [r-^ is a 6m tele- 
scope in the Atacama Desert of Chile, which operated 
at 148, 218, and 277 GHz using three 1024-element ar- 
rays of superconducting bolometers. The maps used in 
this analysis were made over three years of observation 
in the equatorial region during 2008-2010 at 148 GHz, 
and consist of six 3° x 18° patches of sky at a noise level 
of ~ 21 /iK arcmin. In our source mask construction we 
also use maps of the same area made in 2008 at a fre- 
quency of 218 GHz. The maps were calibrated as in |29|. 
We apodize the maps by multiplying them with a mask 
that smoothly increases from zero to unity over 0.1° from 
the edge of the maps. 

Although atmospheric noise is removed in the map- 
making process, we implement an additional filter in 
Fourier space to remove signal at multipoles below ^ = 
500 {£ is the magnitude of the Fourier variable conjugate 
to sky angle). In addition, we remove a stripe for which 
— 100 < idee < 100 along the Fourier axis correspond- 
ing to declination to avoid contamination by scan noise. 
Furthermore, to increase the tSZ signal-to-noise, we ap- 
ply a Wiener filter which downweights scales at which the 
tSZ signal is subdominant. This (non-optimal) filter is 
constructed by dividing the best-fit tSZ power spectrum 
from [isl by the total average power spectrum measured 
in the data maps, i.e. C\^^ /C\^^ . For multipoles above 
£ = 6 X 10^, the tSZ signal is completely dominated by 
detector noise and point sources, and hence we remove 
all power above this multipole in the temperature maps. 
The final Fourier-space filter, shown in Fig. 1, is normal- 
ized such that its maximum value is unity. As it is con- 
structed using the binned power spectrum of the data, it 
is not perfectly smooth; however, we apply the same fil- 
ter consistently to data, simulations, and analytic theory, 
and thus any details of the filter do not bias the interpre- 
tation of our result. After filtering, the edges of the maps 
are cut off to reduce any edge effects that might occur 
upon Fourier transforming despite apodization. Simula- 
tions verify that no additional skewness is introduced by 
edge effects into a trimmed map. 



B. Removing Point Sources 

In order to obtain a skewness signal due only to the tSZ 
effect, any contamination of the signal by point sources 
must be minimized. These objects consist of IR dusty 
galaxies and radio sources. We use two approaches to 
eliminate the point source contribution: template sub- 
traction and masking using the 218 GHz channel. 

In the template subtraction method [30j, which we 
use to remove resolved point sources (mainly bright ra- 
dio sources), sources with a signal-to- noise (S/N) greater 
than five are first identified in a match- filtered map. A 
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FIG. 1. The Wiener filter applied to the ACT temperature 
maps before calculating the unnormalized skewness. This fil- 
ter upweights scales on which the tSZ signal is large compared 
to other sources of anisotropy. 



template with the shape of the ACT beam is then scaled 
to the appropriate peak brightness of each source, and 
this profile is subtracted from the raw data. The pro- 
cess is iterated following the CLEAN algorithm [31] until 
no more sources can be identified. We verify that this 
procedure does not introduce skewness into the maps 
(e.g., through oversubtraction) by checking that similar 
results are obtained using a different procedure for reduc- 
ing source contamination, in which we mask and in-paint 
pixels which contain bright sources with S/N > 5 [32]. 

We take a second step to suppress the lower-flux, unre- 
solved point sources (mainly dusty galaxies) that remain 
undetected by the template subtraction algorithm. At 
218 GHz, dusty galaxies are significantly brighter than 
at 148 GHz and the tSZ effect is negligible. We construct 
a dusty galaxy mask by setting all pixels (which are ap- 
proximately 0.25 arcmin^) to zero that have a tempera- 
ture in the 218 GHz maps larger than a specified cutoff 
value. This cutoff is chosen to be 3.2 times the standard 
deviation of the pixel values in the filtered 148 GHz map 
(3.2cr). This procedure ensures regions with high flux 
from dusty galaxies are masked. We also set to zero all 
pixels for which the temperature is lower than the nega- 
tive of this cutoff, so that the masking procedure does not 
introduce spurious skewness into the lensed CMB distri- 
bution, which is assumed to have zero intrinsic skewness. 
The mask is then applied to the 148 GHz map to re- 
duce the point source contribution. Simulations ([33] for 
IR sources) verify that the masking procedure does not 
introduce spurious skewness into the 148 GHz maps. 

Finally, all pixels more than twelve standard deviations 
(12<j) from the mean are also removed from the 148 GHz 
maps. Due to the ringing around very positive or nega- 
tive pixels caused by the Wiener filter, the surrounding 
eight arcminutes of these points are also masked. This 
additional step slightly increases the S/N of the skew- 
ness measurement by reducing the dependence on large 
outliers, enhances the information content of low mo- 





FIG. 2. Histogram of the pixel temperature values in the 
filtered, masked ACT CMB temperature maps. A Gaussian 
curve is overlaid in red. 



FIG. 3. Likelihood of the skewness measurement described in 
the text (with Gaussian statistics assumed). 



merits by truncating the tail of the pixel probability den- 
sity function, and ensures that any anomalous outlying 
points from possibly mis-subtracted bright radio sources 
do not contribute to the skewness signal. Overall, 14.5% 
of the 148 GHz map is removed by the masking proce- 
dure, though the removed points are random with respect 
to the tSZ field and should not change the signal. 



IV. RESULTS 
A. Evaluating the SkeAvness 

We compute the unnormalized skewness of the filtered 
and processed 148 GHz maps by simply cubing and av- 
eraging the pixel values in real space. The result is 
(T^) = —31 ± 6 /iK^, a 5(7 deviation from the null result 
expected for a signal without any no n- Gaussian compo- 
nents. The skewness of the CMB temperature distribu- 
tion in our filtered, processed maps is visible in the pixel 
value histogram shown in Fig. 2 (along with a Gaus- 
sian curve overlaid for comparison). It is evident that 
the Gaussian CMB has been recovered on the positive 
side by point source masking, with the apparent trunca- 
tion beyond 50 /iK due to the minute probability of such 
temperatures in the Gaussian distribution. The likeli- 
hood corresponding to our measurement of the skewness 
is shown in Fig. 3. 

The "Gaussian statistics assumed" error on the skew- 
ness includes only Gaussian sources of noise. We calcu- 
late this error by using map simulations that consist of 
Gaussian random fields with the same power spectrum as 
that observed in the data, including beam effects. These 
simulations contain Gaussian contributions from IR, SZ, 
and radio sources, the primordial lensed CMB, and de- 
tector noise. This estimate thus does not include the 
error resulting from non-Gaussian corrections, which (af- 
ter source subtraction) are due to the non-Gaussian tSZ 



signal. Though CMB lensing is also a non-Gaussian ef- 
fect it does not contribute to the error on the skewness, 
as the connected part of the six-point function is zero to 
lowest order in the lensing potential, and the connected 
part of the three-point function is also negligibly small 
(see [3^ and references therein). 

We calculate errors that include non- Gaussian correc- 
tions by constructing more realistic simulations. To con- 
struct such simulations, we add maps with simulated tSZ 
signal from ^], which assume as = 0.8, to realizations 
of a Gaussian random field which has a spectrum such 
that the power spectrum of the combined map matches 
that observed in the ACT temperature data. Given the 
simulated sky area, we obtain 39 statistically indepen- 
dent simulated maps, each of size 148 deg . By applying 
an identical procedure to the simulations as to the data, 
measuring the scatter amongst the patches, and appro- 
priately scaling the error to match the 237 deg of un- 
masked sky in the processed map, we obtain a fuh error 
(including non-Gaussian corrections) on the unnormal- 
ized skewness of 14 jaK^. While this error is a robust 
estimate it should be noted that the "error on the er- 
ror" is not insignificant due to the moderate simulated 
volume available. The scatter of skewness values mea- 
sured from each of the simulated maps is consistent with 
a Gaussian distribution. The estimate for the full error 
is coincidentally the same as the standard error, 14 /iK^, 
estimated from the six patches into which the data are 
divided. The full error is used below in deriving cosmo- 
logical constraints from the skewness measurement. 



B. The Origin of the Signal 

Is the skewness dominated by massive clusters with 
large tSZ decrements - as suggested by theoretical con- 
siderations described earlier - or by more numerous, less 
massive clusters? To investigate this question, we mask 
clusters in our data which were found in the 148 GHz 
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FIG. 4. Plot of the skewness signal as a function of the mini- 
mum S/N of the clusters that are masked (this indicates how 
many known clusters are left in the data, unmasked). The 
blue line is calculated using the full cluster candidate catalog 
obtained via matched filtering, while the green line uses a cat- 
alog containing only optically-confirmed clusters 38]. Both 
lines have identical errors, but we only plot them for the green 
line for clarity. Confirmed clusters source approximately two- 
thirds of the signal, which provides strong evidence that it is 
due to the tSZ effect. Note that one expects a positive bias 
of ^ 4 fiK^ for the S/N = 4 point of the blue line due to im- 
purities in the full candidate catalog masking the tail of the 
Gaussian distribution. 



maps using a matched filter as in [36]. All clusters de- 
tected above a threshold significance value are masked; 
we vary this threshold and measure the remaining skew- 
ness in order to determine the origin of the signal. 

Fig. 4 shows a plot of the signal against the cluster 
detection significance cutoff. We include calculations us- 
ing both the full cluster candidate catalog obtained via 
matched filtering as well a catalog containing only clus- 
ters confirmed optically using the methodology of [S^ 
on the SDSS Stripe 82 [381- The SDSS Stripe 82 imag- 
ing data cover ^ 80% of the total map area, and thus 
some skewness signal will necessarily arise from objects 
not accounted for in this catalog. The results for these 
two catalogs agree when masking clusters with S/N > 7, 
but differ slightly when masking lower S/N clusters. This 
effect is likely due to the small shortfall in optical follow- 
up area as well as a small number of false detections (i.e., 
impurity) in the candidate clusters that have not yet been 
optically followed up. 

Using either catalog. Fig. 4 implies that just under 
half of the tSZ skewness is obtained from clusters that 
lie below a 5cr cluster detection significance, while the 
remainder comes from the brightest and most massive 
clusters. The results of [36] suggest that clusters de- 
tected at 5(7 significance are roughly characterized by a 
mass M500 = 5 X 10^^ M^/h^ where M500 is the mass 
enclosed within a radius such that the mean enclosed 
density is 500 times the critical density at the cluster red- 
shift. This value corresponds to a virial mass of roughly 
M = 9 X IO^^Mq/H, which was also found to be the 



FIG. 5. A test for IR source contamination: similar to the 
blue line in Fig. 4, but with a range of values of the cutoff 
used to construct an IR source mask in the 218 GHz band. 
Any cutoff below ^ 3.2a gives similarly negative results and 
thus appears sufficient for point source removal, where a = 
10.3 /jK is the standard deviation of the 148 GHz maps. For 
comparison, the standard deviation of the 218 GHz maps is 
^ 2.2 times larger. The percentages of the map which are 
removed for the masking levels shown, from the least to the 
most strict cut, are 0.7, 2.5, 8.4, 14.5, 23.7, and 36.6%. 



mass detection threshold for high-significance ACT clus- 
ters in [37]. Fig. 4 thus demonstrates that roughly half 
of the tSZ skewness signal is due to massive clusters with 
M > 10^^ M^/h. The theoretical calculations described 
earlier give similar results for the fraction of the signal 
coming from clusters above and below this mass scale, 
which is significantly higher than the characteristic mass 
scale responsible for the tSZ power spectrum signal. 

Finally, the positive value in the full candidate cat- 
alog line shown in Fig. 4 when masking clusters above 
S/N = 4 is consistent with zero. When masking at this 
level (with the candidate catalog which contains some 
impurities), we slightly cut into the negative pixel val- 
ues in the Gaussian component of Fig. 2, leading to a 
small spurious positive skewness. For the points we plot, 
we calculate that this bias is only non-negligible for the 
S/N = 4 cut, where it is ^ 4 /iK^. This bias effectively 
explains the small positive offset seen in Fig. 4. How- 
ever, we discuss positive skewness due to any possible 
residual point source contamination below. Overall, the 
dependence of the measured skewness on cluster masking 
shown in Fig. 4 provides strong evidence that it is caused 
by the tSZ effect. 



C. Testing for Systematic Infrared Source 
Contamination 

Despite our efforts to remove point sources, a small 
residual point source contamination of the signal could 
remain, leading to an underestimate of the amplitude 
of the tSZ skewness. To investigate this systematic er- 



ror source, we vary the level at which point sources are 
masked in the 218 GHz maps (the original level is 3.2 
times the standard deviation of the pixel values in the 
filtered 148 GHz map (3.2<j), as described above). The 
results of this test are shown in Fig. 5, which uses the 
full catalog of cluster candidates as described in Fig. 4, 
since the optically-confirmed catalog does not yet cover 
the entire ACT map. Note that masking at 3.2cr results 
in a skewness measurement which agrees with its appar- 
ent asymptotic limit as the IR contamination is reduced, 
within the expected fluctuations due to masking. While 
a slightly more negative skewness value can be measured 
for some masking levels stricter than the 3.2cr level cho- 
sen in the analysis, fluctuations upon changing the un- 
masked area of the sky are expected, so that it can not 
be rigorously inferred that IR contamination is reduced 
between 3.2(T and 2.65<j. Fig. 5 suggests that masking 
at the 3.2(7 level sufficiently removes any contamination 
by IR sources, and stricter cuts will reduce the map area 
and increase statistical errors unnecessarily. 

However, to further estimate the residual point source 
contamination in the 148 GHz maps, we process simu- 
lations of IR sources from [33| (with source amplitudes 
scaled down by 1.7 to match recent observations, as in 
j39j ) with the same masking procedure as that applied 
to the data (described in §111 B[) . creating a mask in a 
simulated 218 GHz map, and applying it to a simulated 
IR source signal at 148 GHz. We find a residual signal of 
(f^) = 3.9 ± 0.1 /iK^. We treat this result as a bias in 
deriving cosmological constraints from the tSZ skewness 
in the following section. 

We also investigate a linear combination of the 148 and 
218 GHz maps that should have minimal IR source levels, 
namely, an appropriately scaled 218 GHz map subtracted 
from a 148 GHz map. Assuming that the spatial distribu- 
tion of the point sources is not affected by the difference 
in observation frequency between 148 and 218 GHz and a 
single spectral index can be applied to all sources, a sim- 
ple factor of ~ 3.2 [13| relates a point source's signal in 
the two different frequency bands. We find that the ap- 
propriate linear combination (subtracting 1/3.2 times the 
220 GHz map from the 148 GHz map) produces a signal 
in agreement with that resulting from the previously de- 
scribed masking procedure, although the additional noise 
present in the 218 GHz maps slightly reduces the signif- 
icance of the detection. 



V. COSMOLOGICAL INTERPRETATION 

To obtain cosmological information from the measured 
amplitude of the unnormalized skewness, we compare 
our results with two different sets of tSZ simulations 
[33I, [35|- Both sets of simulations are run with ag = 0.8, 
but differ in their treatment of the ICM. The simula- 
tion of [35| is a fully hydrodynamic cosmological simu- 
lation that includes sub-grid prescriptions for feedback 
from active galactic nuclei, star formation, and radiative 



cooling. The simulation also captures non-thermal pres- 
sure support due to turbulence and other effects, which 
significantly alters the ICM pressure profile. The simu- 
lation of [33] is a large dark matter-only A^-body simu- 
lation that is post-processed to include gas according to 
a polytropic prescription. This simulation also accounts 
for non-thermal pressure support (though with a smaller 
amount than [35]), and matches the low-redshift X-ray 
data presented in [17|. 

We perform the same filtering and masking as that 
applied to the data in order to analyze the simulation 
maps. For both simulations, the filtering reduces the sig- 
nal by ~ 95% compared to the unfiltered value. For the 
simulations of [35|, we measure {T^) = — 37 /iK^, with 
negligible errors (the superscript S indicates a simulated 
value). However, this value is complicated by the fact 
that these simulations only include halos below z = 1. 
An analytic estimate for the skewness contribution due 
to halos with z > 1 from Eq. ([T]) gives a 6% correction, 

which yields (T^) = — 39 /jK^. For the simulations of 

|33J . we measure (T^) = — 50 /iK^, with errors also 
negligible for the purposes of cosmological constraints. 

We combine these simulation results with our calcu- 
lated scalings of the skewness and the sixth moment with 
as to construct a likelihood: 
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where (f^) is our measured skewness value and the the- 
oretically expected skewness as a function of as is given 

by 



{fTi^s)={fr{^y 



(4) 



The likelihood in Eq. (|3|) explicitly accounts for the fact 
that a'^Yi^ ^^^ variance of the skewness, depends on as 
— a larger value of as leads to a larger expected vari- 
ance in the tSZ skewness signal. In particular, the vari- 
ance of the tSZ skewness is described by a sixth moment, 
so it scales as a^^ . As determined above, the Gaussian 
and non- Gaussian errors on the skewness are 6 /iK^ and 
Vl4^ — 6^ /iK^ = 12.6 /iK^, respectively. We approxi- 
mate the dependence of the full error on as by assuming 
that only the non-Gaussian component scales with as'-, 
while this is not exact, as some of the Gaussian error 
should also scale with ag, small differences in the size or 
scaling of the Gaussian error component cause negligible 
changes in our constraints on as- 

Finally, although we have argued previously that IR 
source contamination is essentially negligible, we explic- 
itly correct for the residual bias as calculated in the pre- 
vious section. Thus, we replace (T^) = —31 /jK^ with 

(^^)Zrr = -31-3.9 /iK^ in Eq. 1^. (Note that this bias 
correction only shifts the central value derived for as be- 
low by roughly one-fifth of the la confidence interval.) 



Moreover, in order to be as conservative as possible, we 
also model the effect of residual point sources by includ- 
ing an additional IR contamination error (with the same 
value as the residual IR source contamination, 3.9 /iK^) 
in our expression for the variance of the skewness: 

a,\(a8) = 6^ mK^ + 12.62 (^)"^ ^K^ + SJ^ ^K^. (5) 

Using the likelihood in Eq. (|3|), we obtain confidence 
intervals and derive a constraint on as- Our likelihood 
and hence our constraints depend in principle on which 
simulation we use to calculate (T^) , as well as on the 
values we choose for 0^3 and aQ. Using the simulations 
of [351 a nd the scahngs determined above for the profile 
from [ii], we find ag = 0.791^^ (68% C.L.) tom (95% 
C.L.). In Table I, we compare the constraints on as ob- 
tained from the use of different scalings and simulated 
skewness values; the constraints are insensitive to both 
the pressure profile used to derive the scaling laws and 
the choice of simulation used to compute the skewness. 

For comparison, the final release from the Chandra 
Cluster Cosmology Project found as — 0.803±0.0105, as- 
suming ^rn = 0.25 (there is a strong degeneracy between 
as and Qrn for X-ray cluster measurements that probe 
the mass function) [40|. Perhaps more directly com- 
parable, recent studies of the tSZ power spectrum have 
found as = 0.77 ± 0.04 (statistical error only) [13] and 
as = 0.807±0.016 (statistical error and approximately es- 
timated systematic error due to theoretical uncertainty) 
|4H . Our results are also comparable to recent con- 
straints using number counts of SZ-detected clusters from 
ACT and SPT, which found as = 0.851 ± 0.115 (fully 
marginalizing over uncertainties in the mass-SZ flux scal- 
ing relation) [42] and as = 0.807 ± 0.027 (marginalizing 
over uncertainties in an X-ray-based mass-SZ flux scaling 
relation) [43], respectively. Although more than half of 
the tSZ skewness signal that we measure is sourced by de- 
tected clusters (i.e., the same objects used in the number 
counts analyses), our method also utilizes cosmological 
information from clusters that lie below the individual 
detection threshold, which gives it additional statistical 
power. Finally, note that we have fixed all other cosmo- 
logical parameters in this analysis, as as is by far the 
dominant parameter for the tSZ skewness [4^. However, 
marginalizing over other parameters will slightly increase 
our errors. 

To evaluate the theoretical systematic uncertainty in 
the amplitude of the filtered skewness due to unknown 
ICM astrophysics, we test the effect of different gas pre- 
scriptions by analyzing simulations from [35] with all 
forms of feedback, radiative cooling and star formation 
switched off, leading to an adiabatic ICM gas model. For 

these adiabatic simulations we find (T^) = — 56 /iK^ 
(after applying the 6% correction mentioned earlier), 
which for the skewness we measure in our data would 
imply as = 0.77^^;^^ (68% C.L.)to;o5 (95% C.L.). Turn- 
ing off feedback and all sub-grid physics is a rather ex- 
treme case, so the systematic theoretical uncertainty for 



Battaglia 0^3, ae 
Arnaud 0^3, Q^e 
K-S 0^3, Q^6 



r3\S 



Battaglia (f ^) 



n 7Q+0-03 +0.06 

'-'•'^-0.03 -0.06 

n 7Q+0-03 +0.06 

^•'^-0.03 -0.06 



0.79 



+0.03 +0.07 
-0.03 -0.06 



i^3\S 



Sehgal (f ^) 



n 77+0.03 +0.05 

^•' ' -0.02 -0.05 

PI 77+0.02 +0.05 

^•' ' -0.02 -0.05 



0.77 



+0.03 +0.06 
-0.03 -0.05 



TABLE I. Constraints on ag derived from our skewness mea- 
surement using two different simulations and three different 
scalings of the skewness and its variance with as . The top row 
lists the simulations used to calculate the expected skewness 
for (78 = 0.8 [33I, [35] : the left column lists the pressure profiles 
used to calculate the scaling of the skewness and its variance 
with as [il, [13, [2^ . The errors on as shown are the 68% and 
95% confidence levels. 



a typical simulation with some form of feedback should 
be slightly smaller than the statistical error from the 
measurement, though still no n- negligible. This contrasts 
with measurements of ag via the tSZ power spectrum, for 
which the theoretical systematic uncertainty is compara- 
ble to or greater than the statistical uncertainty [13|, |41| . 
As highlighted earlier, this difference can be traced to the 
dependence of the power spectrum amplitude on the ICM 
astrophysics within low-mass, high-redshift clusters. The 
skewness, on the other hand, is dominated by more mas- 
sive, lower-redshift clusters that are less affected by un- 
certain non-gravitational feedback mechanisms and are 
more precisely constrained by observations. Nonethe- 
less, as the statistical uncertainty decreases on future 
measurements of the tSZ skewness, the theoretical sys- 
tematic error will quickly become comparable, and thus 
additional study of the ICM electron pressure profile will 
be very useful. 



VI. CONCLUSIONS 

As the thermal Sunyaev-Zel'dovich field is highly non- 
Gaussian, measurements of non-Gaussian signatures such 
as the skewness can provide cosmological constraints that 
are competitive with power spectrum measurements. We 
have presented a first measurement of the unnormalized 
skewness (T^(n)) in ACT CMB maps filtered for high 
signal to noise. As this is a purely non- Gaussian signa- 
ture, primordial CMB and instrumental noise cannot be 
confused with or bias the signal, unlike measurements 
of the tSZ power spectrum. We measure the skewness 
at 5a significance: (T^(n)) = -31 ± 6 /iK^ (Gaus- 
sian statistics assumed). Including non-Gaussian correc- 
tions increases the error to ± 14 jaK^. Using analytic 
calculations and simulations to translate this measure- 
ment into constraints on cosmological parameters, we 



find as = 0.79+^;^^ (68% C.L.) t^fe (95% C.L.), with 
a slightly smaller but non-negligible systematic error due 
to theoretical uncertainty in the ICM astrophysics. This 
detection represents the first realization of a new, inde- 
pendent method to measure as based on the tSZ skew- 
ness, which has different systematic errors than several 
other common methods. With larger maps and lower 
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noise, tSZ skewness measurements promise significantly 
tighter cosmological constraints in the near future. 
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